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Abstract 

We report an improved measurement of the decay B° — > 7r°7r°, using a data sample of 
535 xW 6 BB pairs collected at the T(45) resonance with the Belle detector at the KEKB 
asymmetric-energy e + e~ collider. The measured branching fraction is B(B° — > 7r°7r°) = 
(1.1 ± 0.3(stat.) ± 0.1(syst.)) x 10 -6 , with a significance of 5.4 standard deviations including sys- 
tematic uncertainties. We also report the partial rate asymmetry: Acp(B° — > 7r°7r°) = 0.44 
^(stat.)i^(syst.). 

PACS numbers: 11.30.Er, 12.15.Hh, 13.25.Hw, 14.40.Nd 
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Measurements of the mixing-induced CP violation parameter sin20i 111, |2j at B factories 
are in good agreement with the Kobayashi-Maskawa (KM) mechanism j3j. To confirm this 
theory, one now has to measure the other two angles of the unitarity triangle, 02 and 03. 
One technique for measuring 2 is to study time-dependent CP asymmetries in B° — > 7r + 7r~ 
decay, where both Belle jij and BaBar (if recently reported the observation of mixing- 
induced CP violation. Belle observed direct CP violation while BaBar found no direct CP 
violation. The extraction of 02, however, is complicated by the presence of both tree and 
penguin amplitudes, each with different weak phases. An isospin analysis of the 7T7T system is 
necessary p, and one essential ingredient is the branching fraction for the decay B° — ■> 7r°7r . 

QCD-based factorization predictions for B(B° — > 7r°7r°) are typically around or below 
1 x 10~ 6 [7], but phenomenological models incorporating large rescattering effects can ac- 
commodate larger values Observation for B° — > 7r 7r° was previously reported by Belle 
with the value of (2.3lo;5(stat.)io;3(syst.)) x 1CT 6 for the branching fraction 0. If such a 
high value persists, an isospin analysis for 2 extraction would become feasible in the near 
future. To complete the program, one would need to measure both the B° and B° decay 
rates, i.e. direct CP violation. 

In this paper we report the improved measurement of the decay B° —> 7r°7r°. We also 
provide a measurement of the direct CP violating asymmetry in this mode. The results are 
based on a 535 xlO 6 BB pairs, collected with the Belle detector at the KEKB asymmetric- 
energy e + e- (3.5 on 8 GeV) collider [[(}. 

The Belle detector is a large-solid-angle magnetic spectrometer that consists of a silicon 
vertex detector (SVD), a 50-layer central drift chamber (CDC), an array of aerogel threshold 
Cerenkov counters (ACC), a barrel-like arrangement of time-of-flight scintillation counters 
(TOF), and an electromagnetic calorimeter comprised of CsI(Tl) crystals (ECL) located 
inside a super-conducting solenoid coil that provides a 1.5 T magnetic field. An iron flux- 
return located outside of the coil is instrumented to detect K\ mesons and to identify muons 
(KLM). The detector is described in detail elsewhere [ill] . Two inner detector configurations 
were used. A 2.0 cm beampipe and a 3-layer silicon vertex detector were used for the first 
sample of 152 xlO 6 BB pairs (Set I), while a 1.5 cm beampipe, a 4-layer silicon detector and 
a small-cell inner drift chamber [l2j were used to record the remaining 383 xl0 6 i?-B pairs 
(Set II). 

Pairs of photons with invariant masses in the range 115 MeV/c 2 < m 77 < 152 MeV/c 2 
are used to form tt° mesons and 7r° mass constraint is implemented; this corresponds to a 
window of ±2.5cr about the nominal 7r° mass, where a denotes the experimental resolution, 
approximately 8 MeV/c 2 . The measured energy of each photon in the laboratory frame 
is required to be greater than 50 MeV in the barrel region, defined as 32° < # 7 < 129°, 
and greater than 100 MeV in the end-cap regions, defined as 17° < 6 1 < 32° and 129° < 
# 7 < 150°, where 9 1 denotes the polar angle of the photon with respect to the positron beam 
line. To further reduce the combinatorial background, 7r° candidates with small decay angles 
(cos#* > 0.95) are rejected, where 9* is the angle between the 7r° boost direction from the 
laboratory frame and one of its 7 daughters in the 7r° rest frame. 

Signal B candidates are formed from pairs of 7r° mesons and are identified by their 
beam energy constrained mass M bc = J E^ 2 am / c 4 — p* B / c 2 and energy difference AE = 
E* B — -Ebeam, where E^ e&ra denotes the run-dependent beam energy, p* B and E* B are the 
reconstructed momentum and energy of the B candidates, all in the e + e~ CM frame. We 
require M bc > 5.2 GeV/c 2 and -0.45 GeV < AE < 0.5 GeV. The signal efficiency is 
estimated using GEANT-based [l3[ Monte Carlo (MC) simulations. The resolution for 
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signal is approximately 3.6 MeV/c 2 in M^ c . The distribution in AE is asymmetric due to 
energy leakage from the CsI(Tl) crystals. 

We consider background from other B decays and from e + e~ — > qq (q = u, d, s, c) 
continuum processes. A large generic MC sample shows that backgrounds from b — > c decays 
are negligible. Among charmless B decays, the only significant background is B ± — > p^ir 
with a missing low momentum tt ± . This background populates the negative AE region, and 
is taken into account in the signal extraction described below. 

The dominant background is due to the continuum processes. We use event topology 
to discriminate signal events from this qq background, and follow the continuum rejection 
technique from our previous publication p. We combine a set of modified Fox- Wolfram 
moments w]| into a Fisher discriminant [15]. A signal/background likelihood is formed, 
based on a Monte Carlo (MC) simulation for signal and events in the M^ c sideband region 
(5.20 GeV/c 2 < M bc < 5.26 GeV/c 2 ) for the qq background, from the product of the 
probability density functions (PDFs) for the Fisher discriminant and that for the cosine 
of the angle between the I?-meson flight direction and the positron beam. The continuum 
suppression is achieved by applying a requirement on a likelihood ratio 1Z sig = £ sig / (Aug + 
£qq), where £ S i g (C q q) is the signal (qq) likelihood. 

Additional discrimination between signal and background can be achieved by using the 
Belle standard algorithm for 6-flavor tagging [Til ], which is also needed for the direct CP 
violation measurement. The flavor tagging procedure yields two outputs: q = ±1, indicating 
the flavor of the other B in the event, and r, which takes values between and 1 and is a 
measure of the confidence that the q determination is correct. Events with a high value of 
r are considered well-tagged and are therefore unlikely to have originated from continuum 
processes. For example, an event that contains a high momentum lepton (r close to unity) 
is more likely to be a BB event so a looser 7Z S i g requirement can be applied. We find that 
there is no strong correlation between r and any of the topological variables used above to 
separate signal from the continuum. 

We divide the data into r > 0.5 and r < 0.5 bins. The continuum background is reduced 
by applying a selection requirement on 7Z S i g for events in each r region of Set I and Set II 
according to the figure of merit (FOM). The FOM is defined as / ^JN^ + A^ g p , where 
and N^ p denote the expected signal, assuming the branching fraction B = 2 x 10~ 6 , 
and background yields obtained from B ± — > p ± 7i° MC and sideband data, respectively. A 
typical requirement suppresses 98% of the continuum background while retaining 45% of 
the signal. 

In our previous publications 0, H3] we reported a special background originating from the 
overlap of a hadronic continuum event and the residual calorimeter energies from an earlier 
QED scattering event. This e + e~ — > 77 background has a signature of back-to-back photons 
in the CM frame, which can combine with two soft photons to form 7r°7r° candidates. The 
fake candidates tend to have very small total CM momenta due to the characteristics of 
the back-to-back photons, resulting in a peaking behavior at M^ c = 5.28 GeV/c 2 . Therefore 
the QED background is quite similar to the B° — > 7r°7r° signal and has to be considered 
in the signal extraction. The best way to identify the e + e~ — > 77 events is to use the 
timing information of the ECL clusters, which are off-time for the QED background but 
only available in the latter 239 fb _1 of data. For data in which the timing information is 
not available, we have to rely on other observables to distinguish candidates synchronized 
with the trigger particles (on-time) from the off-time QED background. Assuming the 
distributions of the off-time QED background are the same for the first 253 fb _1 and latter 
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M hr (GeV/c 2 )' COS0 1 

FIG. 1: (a) The distribution of AE vs. M^ c for the off-time events; (b) the distributions of cos 62 
vs. cos#i for the on-time events (circle) and the off-time events (triangle). 



239 fb _1 , the PDFs of the observables can be obtained using the data in the latter 239 fb _1 
sample, in which trigger timing information can be used to identify the on-time and off-time 
events. 

Figure E^a) shows the M^ c - AE distribution for the off-time 7r°7r° candidates in the 
latter 239 fb _1 sample after passing all analysis requirements. These off-time 7r°7r° events 
are located in the M^ c signal region but have no particular structure in AE. Since the QED 
photons are close to the e + and e~ directions, the angle between the ir° moving direction 
and the z axis can be used to identify the QED background. The angular distributions of 
on-time and off-time candidates in the latter sample are shown in Fig. ^b), where 61 and 9 2 
are the angles of higher and lower momentum 7r°s, respectively. Compared to the on-time 
events, the off-time candidates have a narrow distribution in (cos#i, cos^-plane. The last 
variable that helps distinguish the QED background is missing energy, defined as -E m i SS = 
-^totai " 2 -Ebcam' where E* otai is the total reconstructed energy in the CM frame. Since the 
QED background consists of two overlapping events, the missing energy tends to be larger 
than that of the on-time events. All these five observables are implemented in the fit. 

The signal yields are extracted by applying unbinned five-dimensional maximum likeli- 
hood (ML) fits to the (M bc , AE, cos#i, cos 9 2 , E miss ) distributions of the B and B samples. 
The likelihood is defined as 

C = exp (- £ N. M ) n(E N 'mKJ) C 1 ) 

s,k,j i s,k,j 

where 

Kkj = \^ - <? ■ Acp)PU,kAML x P2 S)fcJ (cos0i, cos^) x P3 sAj (El iss ). (2) 
The direct CP asymmetry is defined as 

N(B -> 7) - N(B -> f) 
ACP ~ N (E^J) + N(B^f)> (3) 
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and s indicates Set I or Set II, k distinguishes events in the r < 0.5 or r > 0.5 bins, j indicates 
the category of signal or background contributions due to the qq continuum, B^ — > p ± n° 
decay and off-time QED events, i is the identifier of the i-th event, Pl Sj fcj(Mb c , AE) are the 
two-dimensional probability density functions (PDFs) in M bc and AE, P2 s ^ k j(cos#i, cos^) 
are the two-dimensional PDFs in cos#i and cos #2) P^s,k,j{E m i SS ) are the one-dimensional 
PDFs in E miss , N Sj ).j is the number of events and q l indicates the B meson flavor: q % = 

+ 1( — 1) for B° and P°. The flavor of the B meson in the B° — > tt tt channel is not self- 
tagged and must be determined from the accompanying B meson. To account for the effect 
of B°-Tr mixing and imperfect tagging, the term Aqp for the signal in Eq. 2 has to be 
replaced by Acpi^- — 2x<f)(l — 2w fc ), where \d = 0.186 ± 0.004 [18j is the time-integrated 
mixing parameter and w k is the wrong-tag fraction that depends on the value of r. The 
wrong-tag fractions are determined using a large sample of self-tagged B° —> D*~tt + , D*~p + 
and D^~l + v events and their charge conjugates |l6| |. 

The M\y C -AE PDFs for the signal and for the B + — > p + 7i° background are taken from 
smoothed two-dimensional histograms obtained from large MC samples. For the signal 
PDF (PI), discrepancies between the peak positions and resolutions in data and MC are 
calibrated using D° — > 7r°7r° and B + — > D°(^ K + n~n°)n + decays. The difference is caused 
by imperfect simulation of the 7r° energy resolution, while the effect of the opening angle 
distributions can be neglected. The invariant mass distribution for the D° is fitted with an 
empirical function for data and MC, and the observed discrepancies in the peak position 
and width are converted to the differences in the peak position and resolution for AE in 
the signal PDF. We require the D° decay products to lie in the same momentum range 
as the 7r°s from B° — > 7r°7r°. To obtain the two-dimensional PDF, Pl(Mb c , AE), for the 
continuum background, we multiply a linear function for AE with the ARGUS function 
for Mb c . The M^-AE PDF for the off-time QED background is modeled as a smoothed 
two-dimensional histogram using the off-time candidates in the latter data set. 

The P2(cos6'i, cos 6*2) and P3(P m i SS ) PDFs are described by two-dimensional and one- 
dimensional smooth histograms, respectively. These PDFs are obtained using the on-time 
and off-time candidates with trigger timing information. Note that the same on-time PDFs 
are used for signals, the p + 7r° background and the continuum. In other words, these three 
components are distinguished based on the Mb c -AP distribution, while the off-time com- 
ponent is identified using all 5 variables. In the fit, the shapes of the signal, off-time and 
B + — > p + 7T° PDFs are fixed and all other fit parameters are allowed to float. We check the 
modeling of the off-time QED background by comparing the results of the five- dimensional 
fit to the data with timing information to the fit result using only the M^ c — AE PDFs after 
removing the off-time candidates. With 3200 events, the obtained off-time yield from the 
5-d fit is 79 d= 14, consistent with 61 off-time candidates. Moreover, the obtained yields for 
signals, the B + — > p + 7r° background and the continuum are also consistent between the two 
fits. 

We perform a 5-d fit to all the data assuming the same off-time QED distributions 
for the former and latter data sets. Figure 121 shows the fit projections. The obtained 
signal yield is 74. 41^^ with a statistical significance (S) of 5.5a, where S is defined as 

S = y / — 2 ln(£ /^7vj) an d £0 an d £>n s denote the maximum likelihoods of the fits without 
and with the signal component, respectively. We vary each calibration constant for the signal 
PDF by ±lcr and obtain systematic errors from the change in the signal yield. Adding these 
errors in quadrature, the systematic error from signal PDF is 

In order to obtain the branching fraction, we divide the signal yield by the reconstruction 
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FIG. 2: Result of the fit described in the text. (Top-left) Mb c projection for events that satisfy 
-0.18 GeV < AE < 0.06 GeV; (top-right) AE projection for events that satisfy 5.27 GeV/c 2 < 
Mbc < 5.29 GeV/c 2 . The bottom plots from left to right are cos6>i, cos #2 and E m i ss projections 
without any M^ c and AE selection. The solid lines indicate the sum of all components, and 
the dotted, dashed, dot-dashed lines and hatched part represent the contributions from signal, 
continuum, off-time and B + — » p + TT° events, respectively. 



efficiency, measured from MC to be 12.8%, and by the number of BB pairs. We consider 
systematic errors in the reconstruction efficiency due to possible differences between data 
and MC. We vary the yields of p^ir and off-time events by ±lcr and obtain the systematic 
error tf.7% and t.2.3%) respectively. We assign a total error of 8% due to 7r° reconstruc- 
tion efficiency, measured by comparing the ratio of the yields of the D — ► K + ir~ and 
D — > K + 7T~7T° decays. The experimental errors on the branching fractions for these de- 
cays |18[ are included in this value. We check the effect of the continuum suppression using 
a control sample of B + —> D°(^ K + 7t~7t )tt + decays; the TZ S i g requirements has a similar 
efficiency for the MC control sample and for signal MC. Comparing the 1Z S requirement on 
the control sample in data and MC, a systematic error of 1.5% is assigned. Finally, we assign 
a systematic error of 1.3% due to the uncertainty in the number of BB pairs, (534.6 ± 7.0) 
x 10 6 , and obtain a branching fraction of 

B(B° -> ttV) = (1.1 ± 0.3(stat.) ± O.l(syst)) x 10" 6 . 

The significance including systematic uncertainties is reduced to 5.4ct. 
This new measurement is lower than our previous published result, (S.Sio^oJ) x 10 — 6 , 
in which the QED background was considered in the systematic uncertainty. We perform 
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FIG. 3: Mb c , AE, cos#i, cos #2 and £W SS distributions with projections of the fit superimposed. 
The distributions are shown separately for events tagged as B° (left) and B° (right). 

a M\y C -AE fit to the data used in the previous analysis. Although the current analysis has 
tighter 1Z s \ g requirements to suppress the continuum background, the obtained value of the 
branching fraction is (2.0ioi) x 10 -6 , consistent with the previous result. We then perform 
a five-dimensional fit to the same data, the central value of the branching fraction drops to 
(1.7±0.4) x 10~ 6 , which is larger than the measurement with 535 million BB pairs. We also 
compare the 2-d and 5-d fit results using the full dataset, the difference of the branching 
fraction is 0.3 x 1CT 6 . Therefore, we conclude that the change in the branching fraction is 
due to a statistical fluctuation and the inclusion of the off-time QED background in the fit, 
where the former has a larger impact. 

The result for A C p is 0.44±g;^±g;g|. Systematic errors are estimated by varying the fitting 
parameters by ±lcr. Including the errors of wrong-tag fraction (u>&) and the time-integrated 
mixing parameter (xd), the total systematic error is li 7 6 7 %. To illustrate this asymmetry, 
we show the results separately for B° and B° tags in Fig. 

In conclusion, we have improved the measurements of B° — > 7r°7r° in a data sample of 
535 million BB pairs. We obtain 74.4^ 9 ; Y signal events with a significance of 5.4 standard 
deviations (a) including systematic uncertainties. The branching fraction is measured to be 
(1.1 ± 0.3(stat.) ± 0.1(syst.)) x 10 -6 . The branching fraction is different from our previous 
result due to the larger data sample that is used and the treatment of the off-time QED 
background. We also report the direct CP asymmetry to be O^Iq^Iq;^, which is consistent 
with our previous result |[|. The branching fraction for B° — ► 7r 7r°, together with the 
measurements of its direct CP violating asymmetry Acp, will allow a model-independent 
extraction of the CKM angle 02 from measurements of the B — > titt system in the near 
future. 
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